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Stability  analyses  of  high-speed  boundary-layer  flow  past  a  5°  half  angle  sharp  cone  with  the 
wall-normal  injection  of  air  through  a  porous  strip  are  performed  using  Navier-Stokes 
solutions  for  the  mean  flow  and  linear  stability  theory.  The  configuration  and  free-stream 
parameters  are  chosen  to  be  similar  to  the  experiments,  which  were  carried  out  at  Caltech’s 
T5  shock  tunnel  to  investigate  the  effect  of  C02  injection  on  laminar-turbulent  transition. 
The  analysis  is  focused  on  pure  aerodynamic  effects  in  the  framework  of  perfect  gas  model. 
It  is  shown  that  the  injection  leads  to  destabilization  of  the  Mack  second  mode  in  the  near¬ 
field  relaxation  region  and  its  stabilization  in  the  far-field  relaxation  region.  To  reduce  the 
destabilization  effect  it  was  suggested  to  decrease  the  injector  surface  slope  or  use  suction¬ 
blowing  of  zero  net  injection.  However,  the  eN  computations  showed  that  these  modifications 
did  not  improve  the  injector  performance  in  the  near-filed  region  in  general.  For  special 
cases  of  low  injection  rates  in  which  the  N-factors  in  the  near  field  region  are  below  the 
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Nomenclature 

=  phase  speed 
=  Mach  number 
=  injection  total  mass  rate 
=  Reynolds  number 
=  pressure 
=  temperature 

=  axial  and  radial  velocity  components 
=  amplification  factor 
=  streamwise  wave  number 
=  azimuthal  wave  number 
=  boundary  layer  thickness 
=  boundary  layer  displacement  thickness 
=  kinematic  viscosity 
=  spatial  growth  rate 
=  angular  frequency 


Superscripts 

*  =  dimensional 


Subscripts 

oo  =  free  stream 
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e  =  upper  boundary-layer  edge 

w  =  wall  surface 


I.  Introduction 

It  is  well  known  that  laminar-turbulent  transition  on  slender  bodies  in  a  hypersonic  flow  at  small  angles  of  attack 
is  caused  by  amplification  of  acoustic  waves  trapped  in  the  boundary  layer  [1,2].  Non-equilibrium  effects  such  as 
molecular  vibration  and  dissociation  can  damp  acoustic  disturbances  [3,4].  Carbon  dioxide  has  been  found  to  be 
well  suited  to  absorb  acoustic  energy  in  an  enthalpy  range  relevant  to  realistic  applications  [5].  These  findings  point 
to  the  potential  beneficial  application  of  adding  carbon  dioxide  into  boundary-layer  flows  in  order  to  delay  transition 
onset.  Experimentally,  this  is  attempted  by  injecting  CO?  into  the  boundary  layer  through  the  wall. 

In  Ref.  [6],  three  injection  schemes  were  studied  both  as  passive  trips  with  no  injection  and  as  active  trips  with 
C02  injection.  The  first  injector  had  four  rows  of  36  orifices.  The  second  injector  was  derived  from  the  former  one 
by  only  keeping  the  fourth  or  most  downstream  row  of  orifices.  The  third  injector  consisted  of  a  micro-porous 
section.  It  was  found  that  the  four-row  injector  tripped  the  boundary  layer  to  fully  turbulent  values  near  the  last  row 
of  orifices.  The  transition  Reynolds  number  with  the  one-row  injector  was  decreased  by  more  than  50%  as  compared 
to  the  smooth  cone.  Therefore,  while  not  as  efficient  as  the  four-row  injector,  one  row  of  holes  still  caused  early 
transition.  A  porous  injector  was  also  tested  and  did  not  result  in  early  transition  when  tested  without  injection. 

Due  to  the  tendency  of  the  discrete  jets  to  cause  transition,  further  studies  were  focused  on  the  porous  injector. 
The  numerical  simulations  [7]  predicted  that  transition  will  occur  immediately  following  the  injection  of  cold  carbon 
dioxide.  The  test  cases  with  air  and  nitrogen  as  the  test  gas  suggested  the  momentum  of  the  injection  plays  a 
dominant  role  in  the  amount  of  amplification  seen  immediately  downstream  due  to  the  inefficient  heating  of  the 
injected  C02.  To  remedy  this,  the  injected  carbon  dioxide  was  pre -heated,  resulting  in  a  reduction  of  amplification  in 
the  post-injection  region.  Despite  this  reduction,  the  stability  analyses  still  predicted  transition  to  occur  earlier  on  the 
cone  as  compared  to  the  case  without  injection.  To  isolate  the  effects  of  non-equilibrium  processes  on  the 
disturbances,  a  second  case  involving  a  similar  cone  with  a  longer  transpiration  interval  was  tested.  For  this  case,  the 
stability  analyses  predicted  a  window  of  carbon  dioxide  injection  resulting  in  a  reduction  of  amplification  and  thus  a 
delay  in  transition.  In  Ref.  [8],  transition  delays  were  documented  in  shots  with  C02  injection.  The  data  showed  a 
general  trend  of  increasing  delay  with  injection  rate,  before  a  sharp  drop-off  at  the  highest  injection  rate.  However, 
the  transition  delay  was  smaller  than  that  observed  in  the  free-stream  with  C02. 

The  foregoing  studies  indicate  that  the  injector  performance  essentially  depends  on  the  balance  between  the 
stabilizing  effect  of  carbon  dioxide  and  the  destabilizing  effect  of  blowing.  The  emphasis  of  this  paper  is  to  study  the 
destabilizing  effects  of  the  injection  process  and  try  to  minimize  them  by  seeking  optimal  distributions  of  the  blow 
rate  and/or  shaping  of  the  injection  region.  The  feasibility  studies  of  different  blowing  configurations  use  the  perfect 
gas  model  as  a  first  step.  Therefore,  the  area  of  most  applicability  from  these  studies  is  near  the  injection  point 
where  aerodynamic  effects  are  most  important  and  the  perfect-gas  model  is  suitable  since  the  injected  gas  is 
relatively  cold  and  not  in  vibrational  or  chemical  non-equilibrium.  Further  downstream,  where  the  injected  flow  is 
heated  and  it  can  effectively  damp  acoustic  instabilities,  the  perfect  gas  model  overestimates  the  disturbance  growth 
rates  and,  thereby,  under-predicts  the  injector  performance.  Once  the  most  promising  injector  configurations  are 
selected  from  this  study,  further  computations  need  to  be  made  of  the  stabilizing  effects  of  C02  downstream  of  the 
injection  region  and  experiments  need  to  be  performed  in  a  high-enthalpy  shock  tunnel  to  confirm  those  predictions. 

The  following  cases  are  considered:  A)  Baseline  configuration  -  a  5°  half-angle  sharp  cone  at  zero  angle  of 
attack  with  normal-wall  blowing  through  a  porous  wall  strip.  This  configuration  is  similar  to  that  tested  at  Caltech’s 
T5  shock  tunnel  [6,8].  B)  Configuration  with  shaping  of  the  injection  region;  C)  Configuration  with  the  normal-wall 
suction  and  blowing  with  zero  net  injection.  For  each  case  the  laminar  basic  flow  is  computed  using  the  in-house 
Navier-Stokes  code  HSFlow,  and  stability  analyses  are  performed  in  the  framework  of  linear  stability  theory  (LST) 
using  the  in-house  LST  code.  These  tools  were  successfully  exploited  for  stability  and  receptivity  studies  of  high¬ 
speed  boundary-layer  flows  [9-12]. 

II.  Baseline  configuration  and  numerical  approach 

As  a  first  step,  we  consider  the  baseline  configuration  shown  in  Figure  1.  A  5°  half-angle  sharp  cone  has  the 
region  128  mm  <  x*  <  169  mm ,  where  air  at  the  wall  temperature  is  uniformly  injected  in  the  normal  wall 
direction.  In  our  computations,  the  injection  total  mass  rate  ranged  from  0  g/s  to  13.5  g/s,  where  the  upper  limit 
approximately  equals  the  double  mass  rate  of  the  incoming  boundary-layer  flow. 
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The  free-stream  parameters  correspond  to  the  free-stream  conditions  of  Run  2540  in  the  GALCIT  T5  shock 
tunnel  [7]:  test  gas  is  air,  Mach  number  M  =5.3,  static  temperature  T*  =  1323.77  K,  static  pressure 

p*x  =  21993  Pa,  density  p*^  =  0.05788  kg/m3,  and  wall  temperature  T*  =  293  K.  The  cone  length  measured 

along  the  cone  axis  is  V  =  1  m.  Hereafter  asterisks  denote  dimensional  quantities.  The  Reynolds  number  based  on 
cone  length  is  Re^  =  4.57  •  106 . 


Figure  1.  A  5°  half-angle  sharp  cone  with  the  injector  -  baseline  configuration. 

Calculations  of  the  basic  flow  were  performed  using  the  in-house  Navier-Stokes  code  HSFlow  for  perfect  gas  of 
Prandtl  number  Pr  =  0.72  ,  specific  heat  ratio  7  =  1.4  and  Sutherland  viscosity-temperature  dependence  with  the 
Sutherland  constant  being  110.4  K.  The  flow  variables  are  made  nondimensional  using  the  free-stream  parameters 
and  the  cone  length  V .  Pressure  is  scaled  using  the  doubled  free-stream  dynamic  pressure  p*JU where  t/*  is 

free-stream  velocity.  The  computational  domain  is  a  rectangle.  Its  bottom  boundary  lies  on  the  cone  surface,  and  the 
upper  boundary  lies  in  the  free  stream  above  the  cone-induced  shock.  The  boundary  conditions  on  the  cone  surface 

are  the  no-slip  condition  and  the  isothermal  condition  T  =  T  .  On  the  inflow  and  upper  boundaries,  the  conditions 

correspond  to  the  undisturbed  free  stream.  On  the  outflow  boundary,  the  unknown  dependent  variables  are 
extrapolated  using  the  linear  approximation. 

The  problem  is  solved  numerically  using  the  implicit  second-order  finite-volume  method  described  in  [9].  The 
two-dimensional  (axisymmetric)  Navier-Stokes  equations  are  approximated  by  a  shock-capturing  scheme  that 
allows  for  modeling  of  flow  non-uniformities  induced  by  the  injector.  The  advection  terms  are  approximated  by  the 
third-order  WENO  scheme  [13]  to  decrease  the  numerical  dissipation.  The  code  algorithm  as  well  as  its 
implementations  and  validations  are  discussed  in  [9].  The  grid  convergence  study  showed  that  the  basic  laminar 
flow  can  be  computed  with  sufficiently  high  accuracy  using  597x649  grid.  This  grid  is  clustered  in  the  direction 
normal  to  the  cone  surface  so  that  the  boundary-layer  region  contains  approximately  50%  of  nodes.  To  resolve  flow 
non-uniformities  in  the  vicinity  of  injector,  the  grid  is  also  clustered  in  the  streamwise  direction  near  the  injection 
boundaries. 

The  mean  flow  profiles  were  used  for  the  local-parallel  stability  analysis.  The  boundary-layer  disturbance  is 
considered  in  the  wave  form  q(y)  exp(iax  +  i/3z  —  iu)t) ,  where  ui  is  angular  frequency  and  a,  (3  are  streamwise 
and  azimuthal  components  of  the  wavenumber  vector,  t  is  time.  The  vector  q  consists  of  the  eigenfunctions  of 
velocity  components  u,v,w  ,  pressure  p  and  temperature  T' .  For  each  x  -station  the  linear  stability  equations  are 
solved  using  a  4th-order  Runge-Kutta  scheme  and  a  Gram- Schmidt  orthonormalization  procedure.  For  the  temporal 
stability  problem,  u}(a,(3,x)  is  a  complex  eigenvalue  while  a,  /3  are  real  parameters.  For  the  spatial  stability 
problem,  a(to,/3,x)  is  a  complex  eigenvalue  while  3,  uj  are  real  parameters.  The  eigenvalues  of  discrete  spectrum 
are  calculated  with  the  help  of  a  shooting/Newton-Raphson  procedure.  The  downstream  amplification  factors 
N(x,lo,/3)  are  computed  using  solutions  of  the  spatial  stability  problem 

X 

N(x,u},(3)=  J  a(u,(3,x)dx ,  (1) 

where  a  =  —a.  is  the  spatial  growth  rate,  and  xQ  is  the  neutral  point  on  the  lower  neutral  branch.  Hereafter  the 

subscript  7’  denotes  the  imaginary  part  of  a  complex  quantity  and  V  denotes  the  real  part.  Stability  computations 
are  performed  with  the  assumption  that  the  second  (bulk)  viscosity  is  zero. 
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III.  Results  for  the  baseline  configuration 

The  mean  flow  fields  of  static  pressure  and  radial  velocity  are  shown  in  Figure  2a, b  for  the  case  of  injection  with 
the  total  mass  rate  m  =  13.5  g/s  (the  nondimensional  mass  flow  rate  is  pwvw  —  (Pwvw)/{P™ —  0.0181).  The 
close  up  views  of  axial  velocity,  static  pressure  and  temperature  fields  in  the  injector  vicinity  are  shown  in  Figure  3, 
where  the  dashed  lines  indicate  the  injection  boundaries.  It  is  seen  that  the  normal-wall  blowing  displaces  the 
boundary  layer  and  induces  compression  waves  emanating  from  the  upstream  boundary  of  the  injector  (red  color  in 
Figure  3b).  A  cold  dead-flow  layer  is  formed  near  the  wall.  Its  thickness  is  of  the  order  of  the  shear-layer  thickness. 
Owing  to  the  mixing  process  this  layer  gets  thinner  and  it  is  ultimately  swallowed  by  the  shear  layer  in  the  mid 
stations  x  «  0.6 .  This  trend  is  clearly  seen  in  Figure  4  where  profiles  of  the  streamwise  velocity  and  temperature 
are  shown  in  various  x  -stations.  Figure  5  shows  the  near-wall  layer  thickness  6(x)  determined  by  the  condition 

U(6)  =  0.99(7  .  After  a  sharp  peak  in  the  injection  region,  <5( x)  decreases  and  slowly  approaches  the  no-blow 

distribution.  These  CFD  data  demonstrate  that  the  injection  strongly  affects  the  near-wall  flow  not  only  near  the 
injector  but  also  in  a  long  post-injection  region  which  is  called  hereafter  as  a  relaxation  region. 

To  get  insight  into  stability  of  the  near- wall  layer  with  injection,  we  start  with  the  temporal  stability  analysis  for 
the  mean  flow  profdes  in  the  station  x  =  0.298  which  lies  in  the  middle  of  the  relaxation  region.  The  flow 
parameters  are  non-dimensionlized  using  their  values  at  the  upper  edge  of  near-wall  layer,  the  length  scale  is  the 

displacement  thickness  <5*  and  the  time  scale  is  6"  /  U * .  In  the  considered  herein  x  -station,  the  local  Mach  number 

is  Me  ss  5  and  the  Reynolds  number  Re^,  =  6*U*  /  v*  w  8074  . 

Using  a  global  search  for  eigenvalues  we  found  seven  unstable  modes  whose  increments  ui.(a)  and  phase 

speeds  cr(a)  =  u>r{a)  /  a  are  plotted  in  Figure  6.  The  computations  were  performed  for  plane  waves  of  /3  —  0. 

Figure  7  shows  the  pressure  eigenfunctions  p(y)  for  modes  0,  1,  2  and  3  with  the  wavenumber  a  corresponding  to 

maximal  increments  lo.  for  each  mode.  These  eigenfunctions  indicate  that  we  are  dealing  with  trapped  acoustic 

waves  propagating  in  a  waveguide  formed  between  the  wall  and  the  relative  sonic  line  y  =  ya  ■  U(y  )  =  cr  —  a(y  ) , 

where  a  is  local  speed  of  sound  [2].  In  the  waveguide  region,  0  <  y  <  ya  ,  where  the  mean  flow  is  supersonic 

relative  to  the  disturbance  phase  velocity,  the  eigenfunction  pn(y)  of  mode  n  oscillates  versus  y  so  that  the 

oscillation  phase  is  approximately  tt  /  4  +  nir ,  n  —  0,1,....  Outside  the  waveguide  (y>ya)  the  eigenfunction 

decays  exponentially.  In  the  Mack  terminology,  mode  0  corresponds  to  the  second  mode,  mode  1  -  to  the  third 
mode  etc.  Because  of  relatively  thick  dead-flow  layer  the  wave  guide  is  essentially  thicker  and  the  number  of 
unstable  acoustic  modes  is  larger  than  in  the  no-blow  case. 
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Figure  2.  Basic  flow  fields  for  the  nondimensional  static  pressure  (a),  and  radial  velocity  (b),  the  injection 
total  mass  rate  m  =  13.5  g/s. 
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Figure  3.  Basic  flow  fields  for  the  axial  velocity,  pressure  and  temperature  in  the  vicinity  of  injector,  the 
injection  total  mass  flow  rate  m  =  13.5  g/s.  The  dashed  lines  show  the  injection  boundaries. 


Figure  4.  Profiles  of  the  streamwise  velocity  U(y)  and  temperature  T(y)  at  various  a; -stations,  the  injection 
total  mass  rate  m  =  13.5  g/s,  the  injection  region  is  shown  by  the  red  line. 
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Figure  5.  The  near-wall  layer  nondimensional  thickness  in  the  no-blow  case  (black  line)  and  in  the  case 
injection  with  m  =  13.5  g/s  (red  line). 


streamwise  wavenumber  a  streamwise  wavenumber  a 

Figure  6.  Increments  tu  (cc)  and  phase  speeds  c  (a)  of  seven  unstable  modes  at  x  =  0.298,  (3  =  0. 
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Figure  7.  The  pressure  eigenfunctions  p  (y)  and  p.(y)  for  modes  0,  1,  2  and  3  at  (cu, a)  corresponding 
maximal  uj.  .  The  red  and  black  lines  show  the  mean- flow  profiles  of  streamwise  velocity  U{y)  and 
temperature  T(y) ,  respectively. 

It  turned  out  that  the  phase  speeds  of  acoustic  instabilities  are  quite  small  compared  with  the  no-blow  case.  As 
shown  in  Figure  8,  the  instabilities  are  observed  in  a  wide  range  of  c  and  their  eigenvalues  are  close  to  the  branches 

of  slow  acoustic  waves  of  the  continuous  spectrum.  The  eigenfunctions  of  such  slow  modes  oscillate  outside  the 
boundary  layer  with  very  weak  damping  (Figure  9).  This  unusual  behavior  may  lead  to  significant  increase  of 
receptivity  to  free-stream  acoustic  disturbances. 

Note  that  maximal  temporal  instability  corresponds  to  mode  2  which  is  the  Mack  fourth  mode  (Figure  6). 
However  the  spatial  stability  analysis  leads  to  the  different  conclusion.  Figure  10  shows  that  the  spatial  growth  rates 
cr(uj)  are  maximal  for  mode  0  corresponding  to  the  Mack  second  mode  -  typical  for  hypersonic  boundary-layer 
flows.  This  switch  over  is  associated  with  abnormal  behavior  of  the  group  velocities  of  the  considered  unstable 
modes. 


Figure  8.  The  phase  speed  trajectories  (red  lines)  of  unstable  modes  in  the  complex  c  plane,  the  branches  of 
continuous  spectrum  are  shown  by  the  black  lines. 
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Mode  0,  x=0.298,  c=0.4 


Figure  9.  Eigenfunctions  p  (y)  and  p  (y)  of  mode  0  (Mack  second  mode)  at  the  phase  speed  c  =  0.4.  The 
red  and  black  lines  show  the  mean  flow  profiles  of  streamwise  velocity  U(y)  and  temperature  T(y) , 
respectively. 


Figure  10.  The  growth  rates  a{w)  and  phase  speeds  c  (w)  of  seven  unstable  modes  at  x  =  0.298,  f3  —  0. 

Figure  1 1  shows  distributions  of  maximal  growth  rates  =  max  x )  for  the  five  unstable  modes  in  the 

relaxation  region.  In  the  all  x  -stations  considered,  the  most  unstable  is  mode  0.  Computations  for  oblique  waves, 
(3^0,  showed  that  their  growth  rates  are  smaller  than  in  the  case  of  f3  —  0 .  These  findings  allow  us  to  restrict 
further  parametric  studies  to  the  Mack  second  mode  of  (3  =  0 ,  which  is  called  hereafter  as  Mack  mode.  Note  that 
the  injection  leads  to  significant  destabilization  of  the  Mack  mode  in  a  relatively  long  region  (compare  the  black  and 
orange  lines  in  Figure  11). 

Distributions  of  the  Mack  mode  growth  rates  (left  plot)  at  various  frequencies  (right  plot)  are  shown  in  Figure  12 
for  the  no-blow  case  (blue  lines)  and  for  the  case  with  injection  of  m  =  13.5  g/s  (black  lines).  The  corresponding 
distributions  of  N-factors  are  shown  in  Figure  13.  As  expected,  the  injection  destabilizes  the  Mack  mode  in  the 
relaxation  region  that  gives  a  local  maximum  of  the  N-factor  envelope  at  x  «  0.3 .  Further  downstream  the  maximal 
growth  rates  remain  higher  than  in  the  no-blow  case.  Flowever  the  width  Ax  of  unstable  regions  is  narrowed  down 
which  leads  to  decreasing  of  the  integral  amplification.  As  a  result  the  N-factor  envelope  is  lower  than  that  of  the  no¬ 
blow  case  in  the  range  0.35  <  x  <  0.75 .  In  the  far-field  region  x  >  0.6  the  mean  flow  is  almost  parallel  and  the 
growth  rates  weakly  depend  on  x .  This  causes  a  steep  rise  of  the  N-factor  envelope.  Eventually  this  envelope 
crosses  that  of  no-blow  case  at  x  ss  0.75 . 
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Figure  11.  Distributions  of  maximal  growth  rates  along  the  relaxation  region  in  the  case  of  injection  with 
m  =  13.5  g/s;  dashed  lines  show  the  injector  boundaries;  the  orange  line  shows  the  baseline  case  without 
injection;  /3  —  0. 


Figure  12.  Growth  rates  and  frequencies  of  the  Mack  mode  for  the  no-blow  (baseline)  case  (blue  lines)  and 
the  case  with  injection  of  m  =  13.5  g/s  (black  lines). 


Figure  13.  N-factors  of  the  Mack  mode  for  the  no-blow  (baseline)  case  (blue  lines),  the  case  of  m  =  13.5  g/s 
(black  lines)  and  the  case  of  m  =  6.75  g/s  (red  lines);  dashed  line  -  N  =  9.2  . 
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The  N  factors  for  the  case  of  m  =  6.75  g/s  are  shown  by  the  red  lines  in  Figure  13.  As  expected,  the  near-field 
maximum  of  their  envelope  is  significantly  smaller  than  in  the  case  of  m  =  13.5  g/s.  However,  the  relaxation 
region  becomes  shorter  and  the  almost  parallel  mean  flow  sets  in  earlier.  The  envelope  of  steep  slope  is  shifted 
upstream  and  the  N-factors  become  larger  than  in  the  no-blow  and  m  =  13.5  g/s  cases  for  x  >  0.5  . 

Thus,  the  considered  injection  leads  to  destabilization  of  the  near-field  relaxation  region,  stabilization  of  the  mid¬ 
field  relaxation  region,  and  destabilization  of  the  far-field  relaxation  region  where  the  basic  flow  is  almost  parallel. 
The  width  and  location  of  these  regions  as  well  as  the  level  of  stabilization/destabilization  effect  depend  on  the 
injected  mass  rate.  To  estimate  the  transition  onset  we  use  hereafter  the  critical  N-factor  N  =  9.2  ,  which  has  been 

derived  in  Ref.  [7]  from  comparison  of  the  eN  computations  with  the  transition  measurements  on  a  5  half-angle  sharp 
cone  tested  in  the  GALCIT  T5  shock  tunnel  without  injection.  For  sufficiently  large  m ,  at  which  the  maximum  of 

N-factor  envelope  N  in  the  near-field  relaxation  region  is  higher  than  the  critical  N  =  Nh  ,  the  transition  onset 
point  xtr  is  predicted  to  be  close  to  the  injector.  Decreasing  of  m  weakly  affects  xtr  until  N  becomes  lower 

than  N .  At  this  moment  the  transition  point  jumps  to  the  a; -station,  which  may  be  even  downstream  from  the 

transition  point  in  the  no-blow  case;  i.e.,  the  premature  tripping  may  be  replaced  by  delay  in  transition.  Further 
decreasing  of  m  to  a  certain  level  can  move  the  transition  point  upstream.  Note  that  the  data  in  Figure  13  (the  first 
and  the  second  regions  of  N-factor  growth)  are  consistent  with  those  of  Ref.  [7]  where  stability  computations  were 
performed  for  CCF  injection  including  real-gas  effects.  Namely,  the  N-factor  envelope  shown  in  Figure  9  of  Ref.  [7] 
for  the  case  of  m  —  13.5  g/s  has  the  same  qualitative  behavior  as  the  envelope  of  black  lines  in  Figure  13. 

IV.  Shaping  of  the  injector 

The  foregoing  stability  analysis  suggests  that  the  injector  performance  may  be  improved  by  shaping  of  the 
injection  region.  It  was  suspected  that  decreasing  of  the  injector  surface  slope  can  partially  compensate  the 
displacement  effect  induced  by  the  blowing  and,  presumably,  reduce  the  instability  growth  in  the  relaxation  region. 
This  speculation  was  examined  by  parametric  stability  studies  of  the  following  configurations. 

The  first  configuration  is  a  conical  shape  schematically  shown  in  Figure  14.  The  rear  part  of  the  original  cone 

(x*  >  169  mm)  is  shifted  along  the  x  -axis  to  A  b" ,  and  it  is  connected  with  the  fore  part  through  a  conical  surface 
(red  line),  which  is  used  for  injection.  The  angle  of  the  new  injector  surface  decreases  from  0  =  5°  to  9  —  A 0  and 

its  x  -length  increases  from  b*  =  41  mm  to  b *  +  A b* .  Stability  computations  were  performed  for  the  injector 
shapes  corresponding  to  A 9  =  1°  and  3°.  The  total  mass  rates  were  m  =  0,  6.75  and  13.5  g/s  for  each  shape. 

The  second  configuration  has  the  injector  of  cylindrical  shape  shown  in  Figure  15.  In  this  case  A 9  =5°,  the 
injector  length  is  41  mm  and  stability  computations  were  performed  at  m  =  0,  1.35,  4.05,  6.75,  9.45  and  13.5  g/s. 


Figure  14.  Conical  shaping  of  the  injection  region. 


Figure  15.  Injector  of  cylindrical  shape. 


A.  Conical  injectors 

Stability  computations  for  the  conical  shapes  of  A 9  =  1°  and  3°  at  zero  blowing  showed  that  the  injector  shapes 
weakly  affect  the  N  factors  of  the  Mack  mode.  As  an  example,  Figure  16  compares  the  distributions  of  N( x)  for  the 

baseline  case  (straight  cone)  with  the  case  of  A 9  =  3°.  These  distributions  are  very  similar  showing  that  the 
shaping  did  not  affect  transition  by  itself. 
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Figure  16.  N-faetors  of  the  Mack  mode  for  the  straight  cone  (black  lines)  and  for  the  shape  of  A9  =  3°  (red 
lines)  at  m  =  0  g/s. 

Now  we  consider  the  baseline  configuration  (A0=O°)  and  the  conical  injectors  of  A9  =  1°  and  3°  at  the 
injection  mass  rate  m  =  13.5  g/s.  Figures  17  and  18  compare  the  basic  flow  fields  for  static  pressure  and  axial 
velocity,  respectively.  As  expected,  the  injector  shaping  reduces  the  pressure  perturbations  induced  by  the  boundary- 
layer  displacement.  However,  the  thickness  of  dead-flow  layer  and,  as  a  consequence,  the  length  of  relaxation  region 
increase  with  increasing  of  A 9  . 
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(c)  x 


Figure  17.  Mean  flow  fields  for  the  nondimensional  static  pressure  at  the  total  mass  flow  rate  m  =  13.5  g/s; 
(a)  -  A9  =  0°  (baseline  cone),  (b)  -  A9  =  1°  ,  (c)  -  A9  =  3° . 
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Figure  19.  N-factors  of  the  Mack  mode  for  the  straight  cone  with  zero  blowing  (blue  lines,  baseline  case),  for 
injectors  of  A 9  =  0°  (black  lines,  cone  with  blowing),  Ad  =  1°  (red  lines)  and  A 9  =  3°  (green  lines)  at 

m  =  13.5  g/s;  dashed  line  -  Nfr  =  9.2  . 

Figure  19  shows  the  Mack  mode  N-factors  computed  for  the  mean  flows  presented  in  Figures  17  and  18.  It  is 
seen  that  the  injector  shaping  produces  a  weak  destabilization  effect  in  the  near-field  relaxation  region 
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0.2  <  x  <  0.4  and  a  noticeable  stabilization  effect  in  the  mid-field  and  far-field  relaxation  regions.  Because  the 
near-filed  maximum  of  the  N-factor  envelope  remains  to  be  higher  than  N  ,  these  injectors  may  cause  premature 
transition. 

B.  Cylindrical  injector 

Figure  20  shows  the  nondimensional  wall-pressure  distributions  P  (x)  of  the  mean  flow  for  different  injection 

rates  m  from  0  to  13.5  g/s.  The  dashed  lines  correspond  to  the  injection  boundaries.  The  baseline  distribution  on  a 
straight  cone  without  injection  is  approximately  uniform  (the  solid  black  line).  The  shaping  without  blowing  (the  red 
line)  leads  to  the  pressure  drop  on  the  injector  surface.  This  behavior  is  consistent  with  the  slender  body  theory  for 

inviscid  flows.  The  blowing  causes  a  sharp  increase  of  P  near  the  upstream  boundary  of  the  injector.  Then  P 
decreases  to  the  level,  which  is  lower  than  in  the  baseline  case.  Further  downstream  the  wall  pressure  is  almost 
constant  in  the  relaxation  region.  The  length  of  P  plateau  increases  while  its  level  decreases  as  the  injection  mass 
flow  rate  increases.  These  mean-flow  changes  affect  the  maximal  (vs.  frequency)  growth  rates  of  the  Mack  mode  as 
shown  in  Figure  21.  The  shaping  itself  (red  line  for  the  case  of  m  =  0)  weakly  affects  the  distribution  of  cr  .  As 

m  increases,  a  sharp  peak  of  (Jn^{x)  is  formed  near  the  upstream  boundary  of  injection  (at  x  «  0.13).  Another 

maximum  is  observed  in  the  relaxation  region.  Because  its  location  is  shifted  downstream  with  m ,  the  injection 
affects  the  N-factor  distributions  in  a  nontrivial  manner  (Figure  22).  The  first  maximum  of  the  N-factor  envelope 
(which  is  located  in  the  near-field  relaxation  region)  increases  with  m  ,  while  the  region  of  steep  envelope  (which  is 
located  in  the  far-field  relaxation  region)  is  shifted  downstream.  Because  the  N-factor  envelopes  are  non-monotonic 

vs.  x* ,  the  eN  method  predicts  different  transition  onsets  depending  on  the  choice  of  critical  N  .  For  N  =  9.2 

(horizontal  dashed  line),  the  most  promising  cases  are  related  to  m  =  6.75  —  9.45  g/s.  More  studies  would  be 
needed  to  account  for  the  non-relaxation  effects  in  the  far  region. 

Figures  23a  and  23b  compare  the  N-factors  for  the  cases  of  Ad  =  0°  (cone  without  shaping)  and  Ad  =  5° 
(cylindrical  injector)  at  rh  =  6.75  g/s  and  13.5  g/s,  respectively.  Although  the  cylindrical  injector  works  better  in 
the  mid-field  and  far-field  relaxation  regions,  its  performance  is  not  improved  in  the  most  critical  near-filed  region 
where  the  N  factor  envelope  has  a  local  maximum. 

Thus,  the  eN  analysis  suggests  that  the  considered  herein  shaping  does  not  prevent  from  early  transition  in  the 
near-field  relaxation  region  at  sufficiently  large  injection  rates  (e.g.,  Figure  23b).  For  relatively  small  m  at  which 
N-factors  in  the  near-field  region  are  below  the  critical  level  (e.g.,  Figure  23a),  the  shaping  produces  a  significant 
stabilization  effect  on  the  boundary  layer  in  the  mid-  and  far-field  regions. 


Figure  20.  The  mean-flow  wall-pressure  distributions  in  the  region  0.1  <  x  <  0.5  at  the  total  mass  rate  from 
0  to  13.5  g/s;  black  line  -  baseline  case  of  straight  cone  without  injection,  dashed  lines  show  the  injection 
boundaries. 
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Figure  21.  Maximal  (vs.  frequency)  growth  rates  of  the  Mack  mode.  Black  line  -  baseline  case  of  straight  cone 
with  zero  injection.  Colored  lines  -  cone-cylinder-cone  configuration  with  different  injection  mass  flow  rates 

from  0  to  mn  =  13.5  g/s.  Dashed  lines  show  the  boundaries  of  injection,  a  is  made  nondimensionai  using 
the  Biasius  length  scale  /  U*  . 


X 


Figure  22.  N-factors  of  the  Mack  mode  for  the  cylindrical  injector  (A  9  =  5°)  at  different  injection  mass  rates 
from  0  g/s  to  mQ  =  13.5  g/s;  dashed  line  -  Nfr  =  9.2  . 
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(a)  (b) 

Figure  23.  N-factors  of  the  Mack  mode  for  the  straight  cone  with  zero  blowing  (blue  lines,  baseline  case),  for 
injectors  of  Ad  =  0°  (black  lines,  cone  with  blowing)  and  cylindrical  injector  with  Ad  =  5°  (red  lines),  (a)  - 
the  injected  mass  rate  is  m  =  6.75  g/s,  (b)  -  m  =  13.5  g/s;  dashed  line  -  N  =  9.2  . 

V.  Suction-blowing  of  zero  mass  injection 

Another  way  to  improve  the  injector  performance  could  be  a  combination  of  normal-wall  suction  and  blowing 
with  zero  net  injection  as  schematically  shown  in  Figure  24.  It  is  expected  that  the  preliminary  suction  of  the 
incoming  boundary-layer  flow  could  partially  compensate  negative  effects  produced  by  the  subsequent  blowing. 


Figure  24  Schematics  of  suction-blowing  of  zero  net  injection. 


This  expectation  is  examined  by  stability  computations  for  the  following  distribution  of  suction-blowing.  In  the 
injection  region  — /*  /  2  <  x*  —  x*  <  T  /  2,  where  x *  =  212.5  mm  is  the  injector  center  and  T  =  41  mm,  the 

mass-flow  rate  q  =  p  v  is  specified  as 

J  •  w  w  1 


q(<f)  =  %  —  sin<A  ^  =  (x 


2t r 


In  this  case  the  suction  mass  flow  rate  is 


x  )2n/l,  —  7T  <  (p  <  7T  . 


the  injection  mass  rate  is 


m 


~%l  /  ^ , 

/  7T  > 


so  that  the  net  mass  flow  rate  is  zero,  in+  +  m_  =  0  .  Computations  were  perfonned  for  m+  =  6.75  g/s  (a  half  of 

maximal  injection  mass  rate  considered  in  the  previous  sections).  The  mean  flow  field  for  static  pressure,  axial 
velocity  and  temperature  are  shown  in  Figure  25.  Figure  26  compares  the  wall  pressure  distributions  for  the  baseline 
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(no-blow)  case,  the  uniform  blowing  of  m  =  6.75  g/s  and  the  suction-blowing  of  m+  =  6.75  g/s.  It  is  seen  that  in 

the  suction-blowing  case  the  relaxation  region  is  shorter  and  the  dead-flow  layer  is  thinner  than  in  the  case  of 
uniform  blowing.  However  the  local  perturbations  in  the  vicinity  of  suction-blowing  system  are  quite  large. 
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Figure  25.  Basic  flow  fields  for  the  nondimensional  static  pressure  (a),  axial  velocity  (b),  and  temperature  (c) 

at  m+  =  6.75  g/s. 


Figure  26.  The  wall  pressure  distributions  for  the  baseline  (no  blow)  case  (black  line),  the  uniform  blowing  of 
m  =  6.75  g/s  (blue  line)  and  the  suction-blowing  of  m  =  6.75  g/s  (red  line). 
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Figure  27  shows  N  factors  corresponding  to  the  cases  in  Figure  26.  As  contrasted  to  the  original  expectation,  the 
suction-blowing  system  destabilizes  the  flow  compared  with  the  uniform  blowing  case.  Thus,  the  preliminary 
suction  does  not  improve  the  injector  performance. 


Figure  27.  N-factors  of  the  Mack  mode  for  the  baseline  (no-blow  case),  the  uniform  blowing  of  m  =  6.75  g/s 
(black  line)  and  the  suction-blowing  of  m+  =  6.75  g/s  (red  line);  dashed  line  -  Nh  =  9.2  . 

VI.  Conclusions 

Stability  analyses  of  high-speed  boundary-layer  flow  past  a  5°  half  angle  sharp  cone  with  the  wall-normal 
injection  of  gas  through  a  porous  strip  were  performed  using  Navier-Stokes  solutions  for  the  mean  flow  and  the 
linear  stability  theory.  The  configuration  and  free-stream  parameters  were  chosen  to  be  similar  to  the  experiments, 
which  were  carried  out  at  Caltech’s  T5  shock  tunnel  to  investigate  the  effect  of  C02  injection  on  laminar-turbulent 
transition. 

It  was  found  that  the  injection  induces  a  cold  dead-flow  layer,  which  is  slowly  swallowed  by  the  shear  layer  in 
the  downstream  relaxation  region.  The  near-wall  flow  behaves  as  a  wave  guide  which  can  support  several  unstable 
modes  of  acoustic  type.  The  spatial  stability  analysis  showed  that  the  most  unstable  mode  corresponds  to  the  Mack 
second  mode.  The  phase  speeds  of  this  instability  are  close  to  the  phase  speeds  of  slow  acoustic  waves  in  the  free- 
stream,  and  its  frequencies  are  several  times  smaller  than  in  the  no  injection  case.  This  may  lead  to  dramatic  increase 
of  receptivity  to  free-stream  noise. 

The  eN  computations  showed  that  the  injection  leads  to  destabilization  of  the  near-field  region,  stabilization  of 
the  mid-field  relaxation  region,  and  destabilization  of  the  far-field  relaxation  region  where  the  basic  flow  is  almost 
parallel.  However,  the  results  obtained  here  are  most  applicable  to  the  near-field  region  where  the  perfect  gas 
assumption  is  valid.  The  width  and  location  of  these  regions  as  well  as  the  level  of  stabilization/destabilization  effect 
essentially  depend  on  the  injected  mass  flow  rate.  These  results  are  consistent  with  the  stability  computations  of 
Wagnild  et  al.  [7],  which  were  carried  out  for  C02  injection  including  the  real-gas  effects. 

In  order  to  improve  the  injector  performance  it  was  suggested  to  decrease  the  injector  surface  slope.  It  was 
suspected  that  this  shaping  could  partially  compensate  the  displacement  effect  induced  by  the  injection  and  thereby 
reduce  the  instability  growth  in  the  near-field  relaxation  region.  However,  the  eN  parametric  computations  predicted 
that  this  shaping  does  not  stabilize  the  near-field  flow  and,  therefore,  it  does  not  prevent  from  early  transition  at 
sufficiently  large  injection  rates.  For  relatively  small  injection  rates,  at  which  the  N-factors  in  the  near-field  region 
are  below  the  critical  level,  the  shaping  produces  a  significant  stabilization  effect  in  the  mid-  and  far-field  regions.  It 
is  expected  that  this  stabilization  will  be  strongly  enhanced  by  the  injection  of  C02. 

It  was  also  suspected,  that  the  injector  performance  could  be  improved  using  a  combination  of  normal-wall 
suction  and  blowing  with  zero  net  injection.  It  was  assumed  that  the  preliminary  suction  of  the  incoming  boundary  - 
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layer  flow  could  partially  compensate  negative  effects  produced  by  the  subsequent  blowing.  However,  the  eN 
computations  showed  that  the  suction-blowing  system  destabilizes  the  flow  in  the  whole  relaxation  region. 

It  should  be  noted  that  the  local  injection  induces  spatial  non-uniformities  of  the  basic  flow,  which  can  alter  the 
instability  growth  rates  via  nonparallel  effects.  Because  our  stability  analysis  does  not  account  for  these  effects,  the 
foregoing  results  should  be  taken  with  great  caution.  The  PSE  analysis  and  direct  numerical  simulations  are  needed 
to  verify  the  eN  predictions  reported  herein.  Note  also  that  the  perfect  gas  model  does  not  account  for  the  CO2 
stabilization  effect  owing  to  absorption  of  acoustic  energy.  This  model  is  suitable  for  the  near-filed  region  where  the 
injected  gas  is  relatively  cold  and  the  absorption  effect  is  weak.  Further  downstream,  where  the  injected  flow  is 
heated,  the  perfect  gas  model  overestimates  the  disturbance  growth  rates  and,  presumably,  under-predicts  the 
injector  performance.  Nevertheless,  the  perfect  gas  modeling  provides  a  good  launching  pad  for  further  stability 
analysis  including  the  real  gas  effects. 
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